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Abstract 

We investigate the vacuum expectation values for the energy-momentum tensor of a mas- 
sive scalar field with general curvature coupling and obeying the Robin boundary condition 
on a spherical shell in the D -\- 1-dimcnsional global monopole background. The expressions 
are derived for the Wightman function, the vacuum expectation values of the field square, 
the vacuum energy density, radial and azimuthal stress components in both regions inside 
and outside the shell. A regularization procedure is carried out by making use of the gener- 
alized Abel-Plana formula for the series over zeros of cylinder functions. This formula allows 
us to extract from the vacuum expectation values the parts due to the global monopole 
gravitational field in the situation without a boundary, and to present the boundary induced 
parts in terms of exponentially convergent integrals, useful, in particular, for numerical cal- 
culations. The asymptotic behavior of the vacuum densities is investigated near the sphere 
surface and at large distances. We show that for small values of the parameter describing 
the solid angle deficit in global monopole geometry the boundary induced vacuum stresses 
are strongly anisotropic. 
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1 Introduction 



The Casimir effect is one of the most interesting manifestations of the nontrivial properties of the 
vacuum state in quantum field theory. Since its first prediction by Casimir in 1948 [1] this effect 
has been investigated for various cases of boundary geometries and different fields (see [2]-[7] and 
references therein). Historically, the investigation of the Casimir effect for a perfectly conducting 
spherical shell was motivated by the Casimir semiclassical model of an electron. It has been 
shown by Boycr [8] that the Casimir energy for the sphere is positive, implying a repulsive force. 
This result was later reconsidered by a number of authors [9]-[12]. More recently new methods 
have been developed for this problem including direct mode summation techniques based on the 
zeta function regularization scheme [13]-[27]. In Refs. [28, 29] the Casimir effect is considered for 
a massless scalar field with the Dirichlet boundary condition on a spherical shell with false/true 
vacuum insidc/outsidc the shell. Recently, the Casimir effect for spherical boundaries on de 
Sitter and anti de Sitter bulks is also investigated [30]-[32]. In the investigations of the Casimir 
effect the calculation of the local densities of the vacuum characteristics is of special interest. 
In particular, these include the vacuum expectation values of the energy-momentum tensor. In 
addition to describing the physical structure of the quantum field at a given point, the energy- 
momentum tensor acts as the source of gravity in the Einstein equations. It therefore plays 
an important role in modeling a self-consistent dynamics involving the gravitational field [33]. 
Investigation of the energy distribution inside a perfectly reflecting spherical shell was made in 
[34] in the case of QED and in [35] for QCD. The distributions of the other components for the 
energy-momentum tensor of the electromagnetic field inside as well as outside the shell, and in 
the region between two concentric spherical shells are investigated in Refs. [36]-[39]. Recently, 
in Ref. [40] the vacuum expectation values of the energy-momentum tensor are evaluated for a 
massive scalar field with general curvature coupling parameter, satisfying the Robin boundary 
condition on spherically symmetric boundaries in Z)-dimensional space. 

The Casimir effect can be viewed as a polarization of vacuum by boundary conditions. An- 
other type of vacuum polarization arises in the case of external gravitational field. In this paper 
we study a situation when both types of sources for the polarization are present. Namely, we 
consider a global monopole spacetime with a concentric spherical shell. Topological defects have 
attracted a great deal of attention because of their relevance to a number of different areas 
ranging from condensed matter to structure formation (for a review see [41]). In the context 
of hot big bang cosmology, the unified theories of the fundamental interactions predict that the 
universe passes through a sequence of phase transitions. These phase transitions might have 
given rise to several kinds of topological defects depending on the nature of the symmetry that 
is broken [42]. If a global SO (3) symmetry of a triplet scalar field is broken, the point like 
defects called global monopoles are believed to be formed. The simplified global monopole was 
introduced by Sokolov and Starobinsky [43]. The gravitational effects of the global monopole 
were studied in Ref. [44], where a solution is presented which describes a global monopole at 
large radial distances. The quantum vacuum effects of the matter fields on the global monopole 
background have been considered in Refs. [45]-[51]. The effects produced by the non-zero tem- 
perature are investigated as well [49]. The zeta function and the Casimir energy for a spherical 
boundary in this background are calculated in [50, 51]. 

In this paper the positive frequency Wightman function, the vacuum expectation values of 
the field square and energy-momentum tensor are investigated for a massive scalar field with 
general curvature coupling parameter ^, satisfying the Robin boundary condition on a spherical 
shell in the D + 1-dimensional spacetime of a point-like global monopole. As special cases they 
include the results for the Diriclet, Neumann, TM, and conformally invariant Hawking boundary 
conditions. To evaluate the corresponding bilinear field products we use the mode sum method 
in combination with the summation formulae from Ref. [52] (see also [53]). These formulae allow 
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(i) to extract from vacuum expectation values the parts due to the global monopolc background 
without boundaries, and (ii) to present the boundary induced parts in terms of exponentially 
convergent integrals involving the modified Bessel functions. We have organized the paper as 
follows. In the next section we consider the vacuum inside a sphere and derive formulae for 
the Wightman function, the expectation values of the field square, energy density and stresses 
on the global monopole background. The behavior of these quantities is investigated near the 
boundary, at the sphere centre and for small values of the parameter describing the solid angle 
deficit in the global monopole geometry. Section 3 is devoted to the vacuum expectation values 
for the region outside a sphere. Section 4 concludes the main results of the paper. 



2 Wightman function, vacuum expectation values of the field 
square and energy-momentum tensor inside a spherical shell 

2.1 Wightman function 

Consider a real scalar field (p with curvature coupling parameter ^ on a -D + 1-dimensional global 
monopole background. In the hyperspherical polar coordinates (r, ??,(?!>) = (r, ^i, ^2, ■ ■ ■ ^n, 0), 
n = D — 2, the corresponding line element has the form 

ds"^ = dt'^ - dr"^ - a^^dnl, (1) 

where dil^ is the line clement on the surface of the unit sphere in D-dimensional Euclidean 
space, the parameter a is smaller than unity and is related to the symmetry breaking energy 
scale in the theory. The solid angle corresponding to Eq. (1) is cr Sd with Sd = 2iT^lyT{D/2) 
being the total area of the surface of the unit sphere in L'-dimensional Euclidean space. This 
leads to the solid angle deficit (1 — (J^)Sd in the spacetime given by line element (1). The field 
equation has the form 

(ViV^ + n? + iR)ip = 0, (2) 

where R is the scalar curvature for the background space-time, m is the mass for the field quanta, 
Vj is the covariant derivative operator associated with the metric corresponding to line element 
(1). The values of the curvature coupling parameter ^ = 0, and = with = {D — 1)/AD 
correspond to the minimal and conformal couplings, respectively. It is convenient for later use 
to write the corresponding metric energy-momentum tensor in the form 



Tik = Vi^VkV ■ 



(3) 



The nonzero components of the Ricci tensor and Ricci scalar for the metric corresponding to 
line element (1) are given by expressions 

I^ = Rl = ... = R^ = n^-^, R = n{n + l)^^, (4) 

where the indices 2,3, correspond to the coordinates 01,02, ■■■ ,4> respectively, and we 

adopt the convention of Birrell and Davies [33] for the curvature tensor. Note that for cr 7^ 1 
the geometry is singular at the origin (point-like monopole), r = 0. 

In this paper we are interested in the vacuum expectation values (VEVs) of the energy- 
momentum tensor (3) on background of the geometry described by (1), assuming that the field 
satisfies the Robin boundary condition 

{Ai + BmWi) cp{x) = (5) 
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on the sphere of radius a, concentric with monopole. Here Ai and Bi are constants, and 
n* = (0, n^, 0, 0) is the unit normal to the sphere, = —1,1 for the interior and exterior regions 
respectively. The imposition of this boundary condition on the quantum field (p{x) leads to the 
modification of the spectrum for the zero-point fluctuations and results in the shift in VEVs for 
physical quantities. In particular, vacuum forces arise acting on constraining boundary. This is 
the familiar Casimir effect. By virtue of Eq. (3) for the VEV of the energy-momentum tensor 
we have 



(0|T,fe(x)|0) = hm did'k{0\<p{xMx')\0) + 



(Ob'(x)|0), (6) 



where |0) is the amplitude for the corresponding vacuum state. Note that the expectation value 
{0\(f{x)(p(x')\0) = G^{x,x') is known as a positive frequency Wightman function. To derive the 
expression for the regularized VEV of the field bilinear product we will use the mode summation 
method. By expanding the field operator over eigenfunctions and using the commutation rules 
one can see that 

{(}\^{xMx')\0) = ^o.{x)ipUx'), (7) 

a 

where {ifaix), (x')} is a complete orthonormal set of positive and negative frequency solutions 
to the field equation with quantum numbers a, satisfying boundary condition (5). 

In the hyperspherical coordinates, for the region inside the sphere the complete set of solutions 
to Eq. (2) with scalar curvature from (4), has the form 



Mx) = (iar-'''^Jui{\r)Y{mk\^,ct>)e-'^\ Z = 0, 1,2, . . . , 
where = (mo = l, mi, . . . , mn), and mi, m2, . . . , rUn are integers such that 

< m„_i < m„_2 < ■■■ <'mi<l, -m„_i < m„ < m„_i. 



(8) 



(9) 



Ju{z) is the Bessel function, and Y{rak]'d, <p) is the surface harmonic of degree I (see [54], Section 
11.2). In Eq. (8) we use the following notation 



i^l 
A 



n 



1 + ^) +{l-a^)nin + l)i^-^D-i) 



.1/2 



(10) 



m^ 



In the following consideration we will assume that uf is non-negative. This corresponds to the 
restriction on the values of the curvature coupling parameter for n > 0, given by the condition 



n 



4(n + l)(c7-2 - 1)' 



(11) 



Note that this condition is satisfied by the most important special cases of the minimal and 
conformal couplings. The coefficients Pa in Eq. (8) can be found from the normalization 
condition 

,2 1 



2a;' 



(12) 



where the integration goes over the region inside the sphere. Substituting eigenfunctions (8), 
and using the relation 



|y(mfe;i?, (P)\^dn = N{mk) 



(13) 
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(the explicit form for N{mk) is given in [54], Section 11.3, and will not be necessary for the 
following considerations in this paper) for the spherical harmonics and the value for the standard 
integral involving the square of the Bessel function, one finds 



2 ^ m,{Xa) 



with the notation 



^'^^^^ " (z2 - ^.2) J2(^) + ^2 j;2(^) • (15) 

From boundary condition (5) on the sphere surface for eigenfunctions (8) one sees that the 
possible values for the frequency have to be solutions to the following equation 

AJ,,{z) + BzJl^{z) = 0, z = Xa, A = A^--B, B = —Bi, (16) 

where = — 1 for the region inside the sphere. It is well known (see, e.g., [54, 55]) that for real 
A, B, and vi > —1 all roots of this equation are simple and real, except the case A/B < —ui 
when there are two purely imaginary zeros. In the following we will assume values oi A/B for 
which all roots are real, A/B > —ui. In terms of the coefficients in (5), it is sufficient to require 



Aia n n 



l + 4(a--l)!i±lc 
n 



(17) 



Note that this condition is satisfied for the most important special cases of Dirichlet scalar, 
Neumann scalar [A/B = 1 — D/2) with the curvature coupling parameter ^ > 0, and for a 
scalar field with the TM type boundary condition [A/B = D/2 — 1). 

Let Xui,k, k = 1,2, ... , be the positive zeros of the function AJ^^(z) + BzJ[,^{z), arranged 
in ascending order. The corresponding eigenfrequencies lo = C0ui,k related to these zeros as 



^i^i,k = \l X\ kl^ ^ Substituting Eq. (8) into Eq. (7) and using the addition formula for 
the spherical harmonics (see [54], Section 11.4), one obtains 



/^-n/2 



(0|^(x)^(a;')|0) = ^^^r__Y,i'^i + n)C';'\cose) (18) 



oo 

X E ^7^f^^^^J.,(A.,.r/a)J.,(A.,,^7a)e-^^'= 



where Cp{x) is the Gegenbauer or ultrasphcrical polynomial of degree p and order and 9 
is the angle between directions (i?, <^) and (??',^'). To sum over k we will use the generalized 
Abel-Plana summation formula [52, 53] 



2Vr,(A,,jk)/(A,,fe) = / f{x)dx+-Kes,=Qf{z)^ 
k=i ^ '^'^ 



'-r dx 



Ux) 



e-^'^7(xe'^^/2) + e'^'^7(xe-'^*/2)l , (19) 



where Iu{x) and Ki,{x) are the modified Bessel functions, and for a given function F{z) we use 
the notation 

F{z) = AF{z) + BzF'{z). (20) 

By taking in this formula i/ = l/2, A = 1,B = 0, asa particular case we receive the Abel- 
Plana formula [56]. Formula (19) was applied previously to a number of Casimir problems 
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for spherically [38, 39, 40] and cylindrically [57, 58] symmetric boundaries on the Minkowski 
background, and for a wedge with and without circular outer boundary [59]. Note that formula 
(19) can be generalized in the case of the existence of purely imaginary zeros for the function 
Jjj {z) by adding the corresponding residue term and taking the principal value of the integral on 
the right (see Ref. [53]). As it has been mentioned earlier, in this paper wc will assume values 
of A/B for which all solutions to Eq. (16) are real. Applying to the sum over k in Eq. (18) 
formula (19), the Wightman function can be presented in the form 

(0|^(x)(^(x')|0) = (0^|(^(a;)^(x')|0^) + {ip{x)^{x'%, (21) 

where 

=0 

and 



(0.|^(x)^(a:')|0.) = ^;^^^^,j:j^,Cr^^^^ J, (.r) J, (./), 

(22) 



1 °° r)7 I 

A'.w'j°.W°)J.-,(^'-7°) 



^Jz'^/a? -m?{t' - t) 



To obtain Eq. (23) we have used the result that the difference of the radicals is nonzero above 
the branch point only. The conditions for the formula (19) to be valid in the case of the sum 
over k in Eq. (18) are satisfied li r + r' + \t — t'\ < 2a. The contribution of the term (22) to the 
VEV does not depend on a, whereas the contribution of the term (23) tends to zero as a ^ oo. 
It follows from here that expression (22) is the Wightman function for the unbounded global 
monopole space with |0m) being the amplitude for the corresponding vacuum. This can be seen 
also by explicit evaluation of the mode sum using the eigenfunctions for the global monopole 
spacetime without a boundary. Note that for a = 1, is equal to / + 1/2 and the sum over / 
in Eq. (22) can be summarized using the Gegenbauer addition theorem for the Bessel function 
[60]. Evaluating the remaining integral involving the Bessel function [61], wc obtain the standard 
expression for the D + 1-dimensional Minkowskian Wightman function. As we have seen, the 
application of the generalized Abel-Plana formula allows us to extract from the bilinear field 
product the contribution due to the unbounded monopole spacetime, and the term (23) can be 
interpreted as the part of the VEV induced by a spherical boundary. In the massless case by 
using the formula for the integral involving the Bessel functions, for the Wightman function (22) 
can be presented as 

,,,,,,,,, 1 ^ (2Z + n)C"/^(cos0)^ {r^ + r'^-{t-t'f + ie\ 
(0.1^(x)^(x')|0.) = ^E ^ (.vtoWU ^ ^-.-V2 ( ^ ) , (24) 

where Qjy, is the Legendre function of second kind, e > for t > t', and e < for t < t'. Note 
that similar expression for the Euclidean scalar Green function is derived in Ref. [46] for = 3 
and in Ref. [48] for an arbitrary D. The corresponding VEVs of the energy- momentum tensor 
are investigated in Refs. [45]-[49]. 

2.2 VEV for the field square 

The VEV of the field square is obtained computing the Wightman function in the coincidence 
limit x' — ^ X. In this limit expression (21) gives a divergent result and some renormalization 
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procedure is needed. As we will see below for < r < a the divergences are contained in the 
first summand on the right of Eq. (21) only. Hence, the renormalization procedure for the local 
characteristics of the vacuum, such as field square and energy-momentum tensor, is the same as 
for the global monopole geometry without a boundary. This procedure is discussed in a number 
of papers (see [45]-[51]) and is a useful illustration for the general renormalization prescription 
on curved backgrounds (see, for instance, [33]). In accordance with this prescription, for the 
renormalization we must subtract from (22) the corresponding De Witt - Schwinger expansion 
involving the terms up to order D. For a massless field the renormalized value of the field square 
has the following general form 

{Om\^\x)\Or,),,r, = [a(''H<7,0 + B(^\a,OHl^r)] , (25) 

where the arbitrary mass scale ^ corresponds to the ambiguity in the renormalization procedure. 
Note that this ambiguity is absent for a spacetime of odd dimension and, hence, B^^^ = for 
even D. Because the dependence of the order of the Legendre function in Eq. (24) on the 
parameter a is not a simple one, it is not possible to obtain closed expression for the coefficients 
A^^^ and B^^\ For small values 1 — approximate expressions are derived in Ref. [48] for 
D = 4: and D = 5. In this paper our main interest are the parts in the VEVs induced by the 
presence of a spherical shell and below we will concentrate on these quantities. 
Using the relation 

. _ r{l + n) 

^' - r(n)/! ' ^^^> 
for the corresponding boundary part in the VEV of the field square we get 

f oo 



where 



is the degeneracy of each angular mode with given I. From the well known properties of the 
modified Bessel functions it follows that in the case of a Dirichlet scalar this quantity is negative 
and for ,f > is a monotonic decreasing function on r. As it has been noted earlier, expression 
(27) is finite for all values < r < a. For a given / and large z the subintegrand behaves as 
^z{rla-\) and, hence, the integral converges when r < a. For large values I, introducing a 
new integration variable y = z/ui in the integral of Eq. (27) and using the uniform asymptotic 
expansions for the modified Bessel functions [60] , it can be seen that to the leading order over I 
the subintegrand multiplied by Di behaves as 

y/l + {zr/ay l + Vl + 2; 

and, hence, the both integral and sum are convergent for r < a and diverge at r = a. For the 
points near the sphere the leading term of the corresponding asymptotic expansion over l/(o— r) 
has the form 

(1 - 2SBo)r (^) 

(-))^ - (4^)P+i)/2(,_,)^-i - (30) 

This term does not depend on the mass and parameter a and is the same as that for a sphere on 
the Minkowski bulk. As the purely gravitational part (25) is finite for r = a, we see that near 
the sphere surface the VEV of the field square is dominated by the boundary induced part, Eq. 
(27). 
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In the limit r — ^ the main contribution into Eq. (27) comes from the 1 = summand with 
the leading term 



r-2(z.o + 1) 



22i/o 



ira 



D-l 



D-l 



f \ 2uQ—n 

a 



f 



dz 



z 



2^0+1 



^^2 _ ^2q2 I^^{z) ' 



0, 



(31) 



where is given by formula (17). As a result the boundary induced VEV (27) is zero at the 
sphere centre for (cr^^ — 1),^ > 0, non-zero constant for (cr^^ — 1)^ = 0, and is infinite for 
(o""^ — 1)^ < 0. As it follows from expressions (25) and (31), for r <C a the total VEV is 
dominated by the purely gravitational part. 

Now we turn to the limit a <C 1 for a fixed value r < a. To satisfy condition (11) we will 
assume that > 0. For > 0, from Eq. (10) one has t'/ ^ 1, and after introducing in Eq. (27) a 
new integration variable y = z/vi, we can replace the modified Bessel function by their uniform 
asymptotic expansions for large values of the order. The integral over y can be estimated by 
using the Laplace method. The main contribution to the sum over I comes from the summand 
with 1 = 0, and to the leading order we receive 



(l-2,5Bo)exp[-2f>ln(a/r)] 
(87ri>)V2r-"5i3C7^-Va2-r2' 



1 



a 



(32) 



For ^ = and cr ^ 1 for the terms with I ^ one has » 1. The corresponding contribution can 
be estimated by the way similar to that in the previous case. This contribution is exponentially 
small. For the summand with Z = to the leading order over a we have vi = nil and, hence. 



Trar^SoO'^ 



ma 



'n/2 



2 2 



e = 0, a<l. 



(33) 



Note that for odd values n the modified Bessel functions in this formula are expressed via the 
elementary functions. In Fig. 1 (left panel) we have plotted the dependence of the boundary 

induced VEV (27) on r/a for the region inside the sphere in the case of a conformally coupled 
(^ = ^d) massless Dirichlet scalar in D = 3 spatial dimensions. The separate curves correspond 
to the values a = 1 (a), a = 0.5 (b), a = 0.2 (c). 





Figure 1: The sphere induced VEV a^^^{ip'^{x))h as a function on r/a in the case of a conformally 
coupled massless D = 3 Dirichlet scalar for the regions inside (left panel) and outside (right 
panel) the sphere. The curves are plotted for cr = 1 (a), cr = 0.5 (b), a = 0.2 (c) 
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2.3 Vacuum energy- momentum tensor inside a sphere 



Substituting the Wightman function (21) with expressions (22) and (23) into Eq. (6), for the 
VEV of the energy-momentum tensor inside a spherical shell one finds 



(0|i;'=|0) = diag (£, -p, -p^, -p^) , 



(34) 



where the vacuum energy density e and the effective pressures in radial, p, and azimuthal, p±, 
directions are functions of the radial coordinate only. Similar to the Wightman function, the 
components of the vacuum energy-momentum tensor can be presented in the form 



qm + ' 



e, P, P±, 



(35) 



where Qm are the corresponding quantities for the monopole geometry when the boundary is 
absent, and quantities qi, are induced by the presence of the spherical shell. As a consequence 
of the continuity equation Vfe(0|r^''|0) = 0, these functions are related by the equation 



r^ + iD-l)(j>,-pu 



0, 



m, b. 



(36) 



For massless fields the VEV of the energy-momentum tensor for the global monopole geom- 
etry without boundaries are investigated in Refs. [45]-[48]. The corresponding renormalized 
components have the structure similar to that given in Eq. (25) for the field square: 



1 



+ ln(/xr) 



(37) 



where the coefficients qm\ Qrn depend only on the parameters a and ^, and Qm^ = in even 
dimensions D. Prom equation (36) the following relations between these coefficients are obtained 



D- 1 



,(1) 



1 



.(2) 



(2) 



D - 1 



(2) 



(38) 



In particular, for conformally coupled fields all coefficients qm ■, qm can be expressed via 

(2) 

£m using the trace anomaly. 

From Eqs. (6), (23), and (27) for the boundary induced parts of the energy-momentum tensor 
components one obtains 



1 oo „ 



r < a, 



(39) 



where for a given function f{y) we have introduced the notations 

'uf 1 + 4^ - 2(mr/?/)2 



Fl^f ifiy)] 



(1 - 40 



f'\y)-^f{y)f'{y) + 



y^ 



1-4^ 



fiy) 



(40) 



^if^ [fiy)] = f'\y) + lf{y)f\y) - ( i + ^^^^1^ ) f\yl i = 4(n + i)C - n (4i) 



4f^Hf{y)] = (4e-i)/''(2/)-^/(y)/'(2/) + 



4^-1 + 



(n + l)y2 



fiy)- (42) 



It can be seen that components (39) satisfy Eq. (36) and are finite for < r < a. 

In the case D = 1 the angular part in line element (1) is absent and we obtain VEVs for the 
one - dimensional segment —a<x<a. In this case (Oni|Tt^|Om)ren = and the boundary part 
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remains only. Due to the gamma function in the denominator of the expression for Di in Eq. 
(28), now the only nonzero coefficients are Dq = Di = 1. Using the standard expressions for 
the functions I±i/2iz) and K±i/2{z) we obtain the formulae for the vacuum energy-momentum 
tensor given in Refs. [40, 62]. Note that in this case the vacuum stresses are uniform and do 
not depend on the parameter ^. For D = 2 and r 7^ the line clement (1) corresponds to the 
flat background spacetime and coincides with the 2D cosmic string geometry. Note that in this 
case n = and from Eq. (10) one has = l/a. 

At the sphere center the main contribution to the boundary induced VEV (39) comes from 
the summands with I = and one has 

, , (2z.o-n)r^'-°-"-V^^^ K,,{z) 
qh(a,r) ^ — ,^ — n — — „ -r / dz , - " / , r — > 0, (43) 

^''^ ^ (2a)2-0+%5D(7-D-lr2(z.o + 1) _ ^2^2 7^^(^) ^ \ > 

where the following notations are introduced 

/o^^^ = (1-40.0, /o^^ = |, /o^-^ = ^^f. (44) 

It follows from Eq. (43) that the boundary induced components are zero at the sphere centre for 
n^((T~^ — 1) > 1, are non-zero constants for n^(c7~^ — 1) = 1, and are infinite otherwise. These 
singularities appear because the geometrical characteristics of global monopole spacetime are 
divergent at the origin. However, note that the corresponding contribution to the total energy 
of the vacuum inside a sphere coming from £{, is finite due to the factor r^~^ in the volume 
element. Comparing expressions (37) and (43) we conclude that near the centre the vacuum 
energy-momentum tensor is dominated by the purely gravitational part. 

Expectation values (39) diverge at the sphere surface, r ^ a. These divergencies are well- 
known in quantum field theory with boundaries and are investigated for various types of bound- 
ary geometries. For the problem under consideration the corresponding asymptotic behavior 
can be found using the uniform asymptotic expansions for the modified Bessel functions, and 
the leading terms are determined by the relations 

Pb^ 2^7r(^+i)/2(a-r)^+i ^ ' ^ ' 

These terms do not depend on mass and parameter a, and are the same as for a spherical shell 
in the Minkowski bulk (see, for instance, [40]). For a conformally coupled scalar the coefficients 
for the leading terms are zero and s,p± ~ (a — r)^^, p ^ {a ~ r)^^^. In general, asymptotic 
series can be developed in powers of the distance from the boundary. The corresponding sub- 
leading coefficients will depend on the mass, Robin coefficient, and parameter a. Due to surface 
divergencies near the sphere surface the total vacuum energy-momentum tensor is dominated 
by the boundary induced parts gj,. 

Now let us consider the VEVs of the energy-momentum tensor in the limit cr <C 1 for a fixed 
r < a. For ^ > by the calculations similar to those given above for the field square, for the 
corresponding boundary parts one receives 

, X _ (l-2JBo)t>3/2exp[-2z>ln(a/r)] (,) 

qi,(a,r) ^ — - /j , (46) 

^ ' (87r)V25Ba«-ir-DV^2372 •'I ' ^ ' 

where D is defined in Eq. (32), and 

/f^ = l-4C, /}"^ = |, fi""^ = D^- (47) 
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Note that in this case the boundary induced vacuum stresses are strongly anisotropic: Pb/pi_b ~ 
C7 <C 1. For a minimally coupled scalar , ^ = 0, the leading term of the asymptotic expansion 
over (T comes from the I = summand in Eq. (39) with vi = ri/2. This term behaves as a^~^ . In 
Fig. 2 (left panel) we have presented the quantities Pb, p±b inside a spherical shell as functions 
on r/a in the case of a conformally coupled massless Z? = 3 Dirichlet scalar on background of 
the global monopole with a = 0.5. Note that for these values of parameters n^((T~^ — 1) < 1, 
and the vacuum energy-momentum tensor components are infinite at the sphere centre. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



Figure 2: The sphere induced VEVs of the energy density, (curves a), radial pressure, 
(curves b) and azimuthal pressure, (curves c) multiplied by a^"*"^, as functions on r/a in the 
case of a conformally coupled massless D = 3 Dirichlet scalar for the regions inside (left panel) 
and outside (right panel) the sphere. The curves are plotted for a = 0.5. 



3 Vacuum densities outside a sphere 

To obtain the VEV for the energy-momentum tensor outside a sphere we consider first the scalar 
vacuum in the layer between two concentric spheres with radii a and b, a < b. The boundary 
conditions on these surfaces we will take in the form 

(^Ai + B,-^^^{x)=0, r = a,b. (48) 

The corresponding eigenfunctions can be obtained from Eq. (8) with the replacement 

J,, (Ar) ^ g,, (Aa, Ar) = J,, (Ar)F,, (Aa) - ^ {Xa)Y,^ (Ar), (49) 

where 1^^ (z) is the Neumann function, and the functions with overbars are defined in accordance 
with Eq. (20). Note that now in the definition of the coefficients in formula (16) we have to 
take = 1. The functions chosen in this way satisfy the boundary condition on the sphere 
r = a. From the boundary condition on r = 6 one obtains that the corresponding eigenmodes 
are solutions to the equation 

C«f (r/, Aa) = J,, (Aa)yjf ) (A6) - J(f) (A6)y,, (Aa) =0, r? = b/a, (50) 

with the notation 

F^'\z) = A,F{z) + B,zF'iz), A, = A,-^, B, = ^, (51) 
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for a given function F{z). In Appendix A we show that for large values rj all roots to Eq. (50) 
are real. As we are interested in the VEVs for the region outside a single sphere, obtained in 
the limit b — > oo, this is the case for our consideration below. Let ji^i^k = Aa , A; = 1, 2, . . . be the 
corresponding positive roots, arranged in ascending order. The coefficients (3a are determined 
from the normalization condition (14), where now the integration goes over the region between 
the spheres, a < r < b. Using the standard formula for integrals involving the product of any 
two cylinder functions one obtains 

^°= 4iv(,„r)L°- ^°-"'-^°'- 

where N{mk) comes from the normalization integral (13) and we use the notation 

= - I [^l + Blin'z^ - ^h] B\z- - ul) I , (53) 

with constants A, B, Ai,, defined in (16) and (51). 

Substituting the eigenfunctions into the mode sum (7) and using the addition formula for the 
spherical harmonics, the expectation value of the field product is presented in the form 

2/ r\—n/2 °° °° 

(0|<^(ar)^(x')|0) = |^g:1^5^(2Z + n)C7,"/'(cos0)5^/i(7.,.)r,f (r?,7.,,fe), (54) 

^ 1=0 k=l 

with the function 



^^i^z^/a'i+nfi{t'-t) 

= / 2 I 2 2 — ^^l"')9^i zr'/a). (55) 



To sum over k we will use the summation formula [40, 52, 53] 



-\"M7.,fc)rr(r?,7.,fc) = / -W^f7^^^^-^-\h{z)\ (56) 
k=l 

TT 



K^u^ ivx) [/i(xe^*/2) + /t(a;e-'^^/2)] dx 
where we have assumed that all zeros for the function (50) are real. In this formula the term 



rv[h{z)\='K ^ Res 

A;,s=l,2 



Hif\z)C-f>{ri,z) 



comes from the residues at the possible zeros z = z^j, of the function Hif\z), s = 1,2, with 
< arg^;] < 7r/2, 

^l-^fc ~ ^l-^fc • Here and below Hu'^\z), s — 1, 2 are the Hankel functions. Note 
that, as we have shown in Appendix A, under the condition (17) the function Ki,{z) has no real 
zeros and, hence, the function Sif\z) has no purely imaginary zeros. In the case of function 
(55) the corresponding conditions for (56) are satisfied ii r + r' + \t — t'\ < 2b. Note that this is 
the case in the coincidence limit for the region under consideration. Applying to the sum over 
k in Eq. (54) formula (56), for the corresponding Wightman function one obtains 

oo 
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zdz Ki'l\r]z) Gui{z,zr/a)G„i{z,zr'/a) ^^^^ 



2 

where we have introduced notation 



^/z'^/a? -m?{t' - t) 



G,{z, y) = h{y)K,{z) - Uz)K,{y), (59) 

with the modified Bessel functions. 

To obtain the Wightman function outside a single sphere on background of the global 
monopole geometry, let us consider the limit 6 — > oo. In this limit the second integral on 

the right of Eq. (58) tends to zero (for large values of the ratio h/a the subintcgrand is pro- 
portional to e~^^^/"), whereas the first integral and the term ry\h{z)\ do not depend on h. It 
follows from here that the quantity 

J p zdz 9uA^,zr/a)g^^{z,zr' /a) jj^2/^2^^2^t,_t-^ \u 
[Jo Vz^ + m^a^ JiM)+nM) '^^^V 

is the Wightman function for the exterior region of a single sphere with radius a. To find 
the boundary induced VEVs we have to subtract the corresponding part for the unbounded 
monopole geometry which, as we saw in previous section, can be presented in the form (22). 
Using the relation 

,.(z, .r/a),.(. .r7a) _ j^^^^/^^j^^^^, ^ _ 1 ^ -^^Hi^\zr /a)Hi^\zr' /a) (61) 
one obtains 

(0|,.(x)^(x')|0) - (0^|^(x),.(x')|0^) = 

, . , -Z) \ H!^Kzr)H(^\zr') + • (62) 

Vz^ + m2 Hlf/iza) ' ' a J 

In this formula we can rotate the integration contour on the right by the angle 7r/2 for s = 1 
and by the angle —tt/2 for s = 2. It can be easily seen that the residue terms coming from the 
poles at the zeros of the functions Si^^ {za) cancel the second summand in the figure braces in 
Eq. (62). Further, the integrals over the segments (0, imo) and (0, —ima) of the imaginary axis 
cancel out and after introducing the Bessel modified functions one obtains 




oo 



i^ixMx')), = J-^^^i±^cr/'(cos^) (63) 



POO 

J ma 



dz J^M) KuAzr/a)K.Mr'/a) ^^^^ 
Kui{z) ^z^-vrP'o^ 



^/z'^/a? ~m?{t' - 1) 



This formula differs from the one for the interior region by the replacements /j/j Ky^ , Ky^ 1^^ . 
For the boundary induced VEV of the field square in the outside region this leads to the formula 



. ^^^^n. r,. (04) 

■Kar-^SoCT^ Jma Ky^{z) ^Jz^-m?a? 
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where Di is defined in accord with Eq. (28). This expression is finite for ah values r > a and 
diverges at the sphere surface. The leading term in the corresponding asymptotic expansion over 
r — a has the form (45) with replacement a — r ^ r — a. In the case of the Dirichlet boundary- 
condition {(p'^{x))h is monotonic increasing negative function on r. 

For the case of a massless scalar the asymptotic behavior of boundary part (64) at large 
distances from the sphere can be obtained by introducing a new integration variable y = zr/a 
and expanding the subintegrand in terms of a/r. The leading contribution for the summand 
with a given I has an order {a/r)^^''^^~^ (assuming that A ^ uiB) and the main contribution 
comes from the I = term. Using the value for the standard integral involving the product of 
the functions [61], the leading term for the asymptotic expansion over a/r can be presented 
in the form 

" - 22^°+i(aa)^-i5^r3(.o + 1)^-^.0 " ^^'^ 

Now comparing this with Eq. (37), we see that for f o > the VEV of the field square at large 
distances from the sphere is dominated by the purely gravitational part. In the limit a the 
leading terms in the asymptotic expansion of Eq. (64) are derived by the way similar to that 
used above for the inside region. For ^ > the corresponding formula is obtained from Eq. (32) 
with replacements — ^ — a? under the square root and a/r ^ r/a under the In function. 
In the case ^ = in Eq. (33) we have to make replacements K^^ , Ki,^ . In the right 

panel of Fig. 1 we have presented the dependence of the boundary induced VEV (64) on r /a for 
the region outside the sphere in the case of a conformally coupled (^ = ^d) massless Dirichlet 
scalar in £) = 3 spatial dimensions. The curves correspond to the values a = 1 (a), a = 0.5 (b), 
a = 0.2 (c). 

As in the interior case, the vacuum energy-momentum tensor is diagonal and the correspond- 
ing components can be presented in the form of a sum of purely gravitational and boundary 
parts, Eq. (35). The parts of these components induced by the presence of a spherical shell are 
given by formulae 



(66) 

where the functions F^f [/(y)] are given by relations (40), (41) and (42). As for the interior 

components, the quantities (66) diverge at the sphere surface, r = a. The leading terms of the 
asymptotic expansions arc determined by same formulae (45) with the replacement a — r — > r — a. 
When D = 1 from Eq. (66) we obtain the expectation values for the one dimensional semi-infinite 
region r > a, given in Refs. [40, 62]. 

For large distances from the sphere, r ^ a, the main contribution into the VEV of the 
energy-momentum tensor comes from the I = summand. The leading terms of the asymptotic 
expansions have the form 



l-OO 

/ dzz^^°-^^Flf[K,,{z)], (67) 
Jo 



where uq is determined by formula (17). Note that the integrals in this formula can be evaluated 
using the value for the integrals involving the product of the functions Ky given in Ref. [61]. 
As we see, for f o > and for large distances from the sphere the vacuum energy-momentum 
tensor is dominated by the purely gravitational part. In Fig. 2 (right panel) we have presented 
the sphere induced quantities e^, p^, p_\_i, for the region outside a spherical shell as functions on 
r/a'm. the case of a conformally coupled massless D = 3 Dirichlet scalar field on background of 
the global monopole with the solid angle deficit a = 0.5. 
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4 Conclusion 



In the present paper we have investigated the Casimir densities induced by a spherical shell 
on background of the D + 1-dimensional global monopole spacetime described by the metric 
(1). The case of a massive scalar field with general curvature coupling parameter and satisfying 
the Robin boundary condition on the sphere is considered. All calculations are made at zero 
temperature and we assume that the boundary conditions are frequency independent. The 
latter means no dispersive effects are taken into account. To obtain the expectation values 
for the energy-momentum tensor we first construct the positive frequency Wightmann function 
(note that the Wightmann function is also important in considerations of the response of a 
particle detector at a given state of motion [33] ) . The application of the generalized Abel-Plana 
formula to the mode sum over zeros of the corresponding combinations of the cylinder functions 
allows us to extract the part due to the global monopole geometry without boundaries and to 
present the sphere induced part in terms of integrals which are exponentially convergent in the 
coincidence limit at any strictly interior or exterior points. The polarization of the scalar vacuum 
by the global monopole without boundaries is investigated in a number of previous papers [45]- 
[49], and our main interest in this paper are the sphere induced parts for the VEVs (Casimir 
densities). The expectation values for the energy- momentum tensor are obtained by applying on 
the corresponding Wightman function a certain second-order differential operator and taking the 
coincidence limit. These quantities diverge as the boundary is approached. Surface divergences 
are well known in quantum field theory with boundaries and are investigated near an arbitrary 
shaped smooth boundary. They lead to divergent global quantities, such as the total energy or 
vacuum forces acting on the sphere and additional renormalization procedure is needed. The 
total Casimir energy for a spherical shell on the global monopole background is investigated in 
Refs. [50, 47] using the zeta function regularization method. 

The boundary induced expectation values of the field square and energy-momentum tensor 
for the region inside a spherical shell are given by formulae (27) and (39) respectively. These 
expressions are finite at interior points and diverge on the sphere surface. The leading term in 
the corresponding asymptotic expansions is the same as for a sphere in the Minkowski bulk and 
is zero for a conformally coupled scalar. The coefficients for the subleading asymptotic terms will 
depend on the boundary curvature, Robin coefficient, mass and the parameter a describing the 
solid angle deficit. In section 3 to deal with discrete modes, first we consider the scalar vacuum 
in the spherical layer between two surfaces. The quantities characterizing the vacuum outside 
a single sphere are obtained from this case in the limit when the radius of the outer sphere 
tends to infinity. Subtracting the parts corresponding to the space without boundaries, for the 
boundary induced parts of the Wightman function and vacuum densities we derive formulae 
(63), (64) and (66). They differ from the ones for the interior region by the replacements 
lu^ , Kin ~^ -^t^i • Near the sphere centre and at large distances from the sphere the VEVs 

for the field square and energy-momentum tensor are dominated by the purely gravitational 
parts corresponding to the monopole geometry without boundaries. For the points near the 
sphere surface the boundary induced parts in the VEVs dominate. We have also investigated 
the VEVs in the limit of small values for the parameter cr, describing the solid angle deficit, 
£7^1. In this limit, corresponding to strong gravitational fields, for ^ > the boundary induced 
parts behave as (v?^)^ ~ cr^/^"-^ exp(— 7/(7) and £{, ~ p_Lb ~ Pb/c ~ cr"'^'"''"''^ 6xp(— 7/(7), with 
7 = 2-\J n{n + 1).^| ln(a/r)|, and the boundary induced stresses are strongly anisotropic. For a 
minimally coupled scalar (^ = 0) and cr <C 1 we have shown that (</3^)6 ~ % ~ ■ <1 = £-P±'P- 
As an illustration of general results, in Fig. 1 and Fig. 2 we have plotted the sphere induced 
VEVs of the field square and energy- momentum tensor components as functions on r/a for a 
conformally coupled Dirichlet scalar in I? = 3 spatial dimensions. Note that in this case the 
leading terms in the asymptotic expansions of the vacuum energy-momentum tensor components 
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near the boundary vanish. The VEV of the field square is negative everywhere and the energy 
density is negative/positive inside/outside the sphere. 
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A On the zeros of the function C"^(?7, z) 

Here we will first show that for real values of the coefficients A, B, A^, Bb in (16) and (51) all 
roots of the equation 

Cfiv, z) = Mz)Y(''^ ivz) - i;(z)4^) ivz) = 0, r? > 1 (68) 

are real or purely imaginary. To see this note that from the relation (C"''(ry, z))* = C^'^{r], z*) it 
follows that if z is a root to (68) when z* is also a root. Further, from the Bessel equation for 
the function gu{z, zx) defined in (49) (considered as a function on x) the following orthogonality 
relation can be derived 



dxxg,,{zi,zix)g,y{z2,Z2x) = 0, zl ^ zl, (69) 



'1 

for two roots zi and Z2 of Eq. (68) . Assuming / z\ and taking in this relation Z2 = z\, we 
obtain a contradiction: 

n 

dx x\gy{zi, zix)\ =0, (70) 

'1 

and, hence, z'f' = z\. This means that the zeros of the function C^^{r],z) are real or purely 
imaginary. 

In this paper we are interested in the limit r] ^ oo. Let us show that for large values rj all 
zeros are real. For the purely imaginary z, z = ye'^*/^ with real y, the function C^^{rj, z) can be 
expressed via the Bessel modified functions: 

C'(^,ye-/2) = - \K,{y)ii'\vy)-Uy)Ki'\w)\ ■ (7i) 

TT L J 

For large rj the contribution of the second term in the square brackets is exponentially suppressed 
and one has 

(r?,2/e-/2) « -i^,(y)lW(TO), V » 1- (72) 

TT 

From the asymptotic expansion of the function Iv{z) for large values of the argument it can be 
seen that for large rj the function I^\riy) has no zeros. Further, using the recurrent relations 
for the Bessel modified functions and expressions (16) for the coefficients A and B with = 1, 
we can see that 

(73) 



a 



^ + i^-^]Ku{y) + yK,_^{y) 



Under the condition (17) with v > the expression in the square brackets of this equation is 
always positive and the function K^{y) has no real zeros. As a result we conclude that for large 
rj the function C^^{rj,ye^^/'^) has no real zeros with respect to y. Therefore, for sufficiently large 
rj all zeros of the function C^{rj, z) are real. 
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